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^\. , We derive upper bounds for spin velocity in half-integer-spin Heisenberg antiferromagnetic chains. 

O^ ' We relate these upper bounds to the instability of the gapless state, which is observed in frustrated 

^\ ' systems. 
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^^ ' During the last decade low dimensional strongly correlated electronic systems have become a subject of intensive 

. I experimental and theoretical research (see M and references therein) . An important class of such systems are one- 
dimensional electronic chains. It is believed, that systems with an odd number of electrons per unit cell (and, therefore, 

'~~ i. half-integer spin per unit cell) can be in a gapless phase. Such a phase should have spin-1/2 fermionic excitations 

1^ ' - spinous, with a linear dispersion law for small values of the wave vector k, E = ck. In the integrable case of the 

5_i ' spin-1/2 nearest-neighbor Heisenberg model existence of such excitations was established by Faddeev and Takhtajan 

"^ ' in the paper [Ej . In the nonintegrable case no proof of the above picture exists, but its correctness seems plausible. 

1 t I In this paper we consider one-dimensional antiferromagnetic Heisenberg models with a half-integer spin per unit 

2 ' cell. Taking for granted the picture of gapless excitations with a linear dispersion law for small k, E = ck, we derive 
r^ [ an upper bound on the spin velocity c. It is then shown, that the condition c > puts a limit on the stability of the 

' , ■ gapless phase in the case, when the system is frustrated. Finally, it is shown on examples, that the derived upper 

-H ' bound is rather good, and in the spin-1/2 case works much better, than the linear spin-wave approximation. 

Q Consider a one-dimensional half-integer-spin Heisenberg model, which is defined on a ring. The Hamiltonian of the 

O model reads 

^; H = ^j,(s„-s„+z). (1) 

>. ! n,l 

QQ Below we argue that the following upper bound holds for the spin velocity of spinous c. 

O ; c < - ^(4/3);Vi7r(S„ • S„+;) (2) 

a^: 

0^ I Here (S„ • S„+;) denotes the spin-spin correlation function. Since we use periodic boundary conditions, this function 
does not depend on the index n. We note, that for the nearest-neighbor model the nonequality (y) relates the spin 
velocity to the energy per unit cell. 
G ' Here we describe a sketch of the derivation for the case of the nearest-neighbor 5" = 1/2 Heisenberg model. The 

i-rt ■ full derivation is given at the end of the paper. 
^ ' 1. Consider a periodic ring of 2N + 1 sites. The full spin of the ground state | ^q ) is 1/2, which corresponds to 

^ having one spinon in the fc = state. Let us denote the energy of the ring as EJ^ ' . 

2. We apply Lieb-Schultz-Mattis transformation pj to this state, obtaining a new quantum state | i/'i )• We then 
check, that this state has a wave vector k = 7r(l-|-l/(2Af+l)). Using this state as a trial state, and evaluating its energy, 
^^ ' we show that in the thermodynamic limit the energy difference between the lowest energy state at /c = 7r(l-|-l/(2iV-t-l)) 
and the ground state of the system is less or equal to 
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Ei-E„< -{i/3)j7,'{S„ ■ S„+,)/(2JV + 1). (3) 

3. Since the lowest energy spectrum E{k) is it periodic, as was proven in B, one can evaluate the spin velocity as 
[E{Tr + Ak) - E{7t)] /Afc, where Ak is 7r/{2N + 1). Then, using the fact, that E{tt) > Eg {Eg is the ground state 
energy) and (m, one obtains (0). 

Let us now study the implications of the inequality (0). For the nearest-neighbor model the inequality (0) gives 
c < —{4:/3)TTEper.site — 0.59l7rJ. The exact value of c given by the Bethe Ansatz is O.SttJ. Let us note that the 
linear spin-wave theory gives c ~ O.BlSvrJ. This shows, that, for spin-1/2 models, the wave function of spinous is 
much closer to the modes, described above, than to the linear spin waves. We also note, that the inequality (pi) may 



be improved, if one finds a way to modify the LSM transformation in such a way that it does not break the rotation 
group symmetry. Studies in this direction will be reported elsewhere B 

Consider now a frustrated S = 1/2 chain with antiferromagnetic exchanges Ji,J2- Then since the spin- velocity 
should be positive, the gapless state is necessarily unstable for Ji(S„ • S„+i) + 4J2(S„ • S„+i) > 0, i.e., for 

J2/J1 > -0.25(S„ • S„+i)/(S„ • S„+2). (4) 

Using correlation functions ratio from the unfrustrated system we get a rough estimate of the critical value of J2/J1 
at 0.5. Numerical calculations [^ put this value at about 0.4. At the transition point the instability leads to an 
appearance of a dimerized phase with a gap. Let us now consider the largc-S limit. In the limit of large S the 
groundstate is the Neel state and then we get the condition J2/J1 > 0.25. Let us note that J2 = 1/4 is the point, 
where the classical Neel state becomes unstable with respect to spiral order, so the classical limit is recovered exactly. 

Derivation 

Here we use notations similar to Q. We first give a derivation for the case of the nearest-neighbor S = 1/2 model, 
and then show how it is extended to the general Hamiltonian (nf). 

l.We assume that the system is in the gapless antiferromagnetic phase, and that the ground state of a ring having 
2N+1 sites has the full spin S^^i^ = 1/2. We denote the ground state of the ring as | i/'o )• We choose the coordinate 
system in such a way, that the z-component of the full spin S? 11 is equal to 1/2. 

'full 
Any Sr 1, = 1/2 configuration | a) is characterized by A^ + 1 numbers a;i < 12 < X3... showing the positions of 

fictitious particles, each of them increasing S^ on a given site by one. The vacuum of the system is assumed to have 

all spins down. Since the system is periodic, the x's are defined only mod 2N + 1. Consider the state 



2. Let us write the state | i/'o ) ^^ ^ linear combination of S? 11 = 1/2 spin configurations | "00 ) — 'J2a^c I ^) 



I V-J^) >= ^ e5*^("^+"^+-)A, I a). (5) 

Adding 2A^ + 1 to any of the x's does not change the \ ip{ >. Therefore, the mentioned freedom of defining the a;'s 
is satisfied. If the translation operator operator T : n ^ n + 1 acts on the state \ i/jl >, then all a;'s are incremented 

by 1. Therefore, an additional phase factor e"[^+^/(^^+^)l is acquired, and the state | -01 > has a wave vector 
7r[l + l/(2Ar+l)]. 

Let us evaluate the expectation value of energy {ipi \ H \ ^p[ ) . The Ising part of the Heisenberg Hamiltonian 
does not flip spins. In contrast the exchange part may change the position of one spin-up fictitious particle, which 
amounts to a phase factor e ^"+1 , depending on a direction of the move. This explains the result of a formal 
calculation, which gives: 



-E-/ 



(0W I H I 0W) = (6) 

[1/2 (e^(v.w I s^s-^, 1 4^)) +e-^(^r^ I st^.s- \ ^r^>) + (4''^ I ^n+i^^ I v^rv . 



Taking a thermodynamic limit of (o) one finds 



(4"^ I H I 0^)) - (4"^ I H I 4"^) - (7) 



J, 



',(^) I Q+ Q- I ,/,W\ ^.^^rrT/,/,W I Q+Q- I ,/,W\ 



{2N+l)-ie~^^{^^"> I S:^,S- I ^^'")+e^^(^l,"> \ S^S'^, \ V^^ 



^ ((^""^ I S^S-^, I 4^') + (4^' I 5+^,5- I V^r^)) + 0(l/(2A^ + 1)^) 



2N- 



The operator S^S^_^_i — S^_^_iS~ is proportional to the current operator of the z-component of the spin [Q . Since the 
ring has reflection symmetry, one can always choose a ground state in such a way, that the expectation value of this 
operator is exactly zero. Then one has 

(4"^ I H I V4"^) - (4"^ I H I V.r^) - (8) 

(^r^ I sis-^, 1 4^)) + (0^) 1 5+^,5- 1 vr^)) + o{i/{2N + If) = 

= -(4/3) j7r2(S„ . S„+i)/(2A + 1) + 0(1/(2A^ + 1)^). 



2A + 1 



3. Now we note, that, as was proven in [Q, the spectrum of the lowest energy excitations E{k) is 7r-periodic in the 
thermodynamic hmit. Therefore, one can find the spin-velocity using the formula c = [£{11 + Afc) — -^(Tr)] /Afc, where 
Afc is Afc = tt/{2N + 1). The energy of the lowest energy state at fc = tt + Afc is less or equal to the expectation value 
of energy for the trial state | ■0i ) ■ The energy of the lowest energy state at fc = tt is more or equal to the energy 
of the ground state (equality is reached in the thermodynamic limit). Based on these considerations, taking (||) and 
dividing it by Afc we obtain in the thermodynamic limit iV ^ 00 an upper bound for the spin velocity c 

c < -(4/3)j7r(S„ • S„+i) (9) 

Now for the general case O) , repeating calculations above and introducing a sum over the index I in the Hamiltonian, 
one arrives at the general inequality (0). One can also repeat the calculations for the case of a general half- integer 
spin s, again obtaining (ph. Let us note that for large S the upper bound (0) scales as S*^, whether the linear-spin 
wave approximation, which becomes exact in large-S limit, gives the spin velocity, which scales like S, in particular 
for the nearest-neighbor model one has c — 2SJ. For the S — 3/2 nearest-neighbor model numerics give c = 3. 87 J, 
linear spin waves give c — 3 J, and the upper bound (y) is 11. 8 J (we use the exact value of the ground state energy 
from 01). This shows, that the upper bound derived here is interesting only for S — 1/2 chains. 

To conclude, we have derived an upper bound for the spin velocity of half-integer-spin Heisenberg antiferromagnetic 
chains and have shown its relation to the instability of the gapless state in the presence of frustration. 
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